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Abstract Radio-quiet 7-ray pulsars like Geminga may 
account for a number of the unidentified EGRET sources 
in the Galaxy. The number of Geminga-like pulsars is 
very sensitive to the geometry of both the 7-ray and ra- 
dio beams. Recent studies of the shape and polarization 
of pulse profiles of young radio pulsars have provided ev- 
idence that their radio emission originates in wide cone 
beams at altitudes that are a significant fraction (1 -10%) 
of their light cylinder radius. Such wide radio emission 
beams will be visible at a much larger range of observer 
angles than the narrow core components thought to orig- 
inate at lower altitude. Using 3D geometrical modeling 
that includes relativistic effects from pulsar rotation, we 
study the visibility of such radio cone beams as well as 
that of the 7-ray beams predicted by slot gap and outer 
gap models. From the results of this study, one can ob- 
tain revised predictions for the fraction of Geminga-like, 
radio quiet pulsars present in the 7-ray pulsar popula- 
tion. 

Keywords pulsars • 7-ray sources • pulsar populations • 
non-thermal radiation 
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1 Introduction 

Rotation-powered pulsars and their winds are presently 
the only known Galactic gamma-ray sources emitting at 
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energies above 100 MeV. The EGRET detector on the 
Compton Gamma-Ray Observatory (CGRO) detected 
six 7-ray pulsars with high confidence [44]. five of which 
were known radio pulsars. The sixth, Geminga, is radio 
quiet (or at least an extremely weak radio source [30]). In 
addition to detecting many sources with known counter- 
parts at other wavelengths, EGRET detected some 110 
sources [5] with no known counterparts, about a third of 
which seem to be of Galactic origin. It is quite plausi- 
ble that a sizeable fraction of these unidentified EGRET 
sources are radio-quiet or radio-weak 7-ray pulsars, i.e. 
Gemingas. There are two possible reasons that may cause 
a 7-ray pulsar to be radio-quiet. The 7-ray emission may 
have a wider beam than the radio beam or be oriented 
in a different direction, or the radio emission may be ei- 
ther absent or too weak to be detectable by current radio 
telescopes. Geminga-like pulsars are therefore very likely 
candidates for some of the unidentified EGRET sources 
in the Galactic plane. 

It has proven difficult to discover Geminga-like pul- 
sars since EGRET typically did not collect enough pho- 
tons in a source to allow pulsation searches without an 
ephemeris known from other wavelengths. Geminga it- 
self, the second brightest EGRET source, was only iden- 
tified as a 7-ray pulsar after its period was discovered by 
ROSAT in the X-ray band [25], although a direct detec- 
tion in EGRET data was later found to be feasible [4]. It 
is therefore useful to estimate the number of Gemingas 
expected through modeling of the 7-ray pulsar popula- 
tion assuming different radio and 7-ray emission models. 
A number of studies of this type have been done ( [17] 
[j"8] for the polar cap model and [9] and [27] for the outer 
gap model) . The results of these studies have shown that 
the Geminga fraction, or the number of radio-quiet 7-ray 
pulsars relative to the total detectable number of 7-ray 
pulsars, strongly depends on the 7-ray emission model. 

There are two main types of 7-ray emission models 
for which population studies have been carried out. Polar 
cap models assume that particle acceleration takes place 
near the neutron star magnetic poles and that 7-ray 
emission results from cascades initiated by magnetic pair 



2 



Alice K. Harding et al. 



production in strong magnetic fields [13] [12] • Outer gap 
models assume that the particle acceleration takes place 
in vacuum gaps that form in the outer magnetosphere[8] 
and that 7 rays result from cascades initiated by pair pro- 
duction of 7 rays and soft X-ray photons from the neu- 
tron star polar cap [7] [39]. The population studies that 
assume these two types of models show that the Geminga 
fraction in outer gap models is much higher than in polar 
cap models. However, most of these studies [17] [9] have 
assumed very simplified geometry for the radio and 7- 
ray emission beams, taking constant solid angles for the 
radio and 7-ray beams, or assuming random and inde- 
pendent relative orientations for these beams. Gonthicr 
et al. [18] [20] have more recently explored more complex 
and realistic emission models for the radio and polar- 
cap/slot-gap 7-ray beams. Jiang et al.[57] have stud- 
ied the solid angle dependence of outer gap 7-ray emis- 
sion beams in their population model, but did not con- 
sider the relative orientations of the radio and 7-ray 
beams. Furthermore, there has never been a single study 
that compares polar-cap/slot-gap models with outer gap 
models on the same footing, using an identical set of as- 
sumptions for the neutron star population, its evolution 
and radio emission geometry, taking into account the rel- 
ative orientations of the radio and 7-ray beams. 

The description of the radio emission in these stud- 
ies has necessarily been empirical, since a physical model 
for the coherent radio emission does not exist. The em- 
pirical models that have been developed from study of 
radio pulse profile morphology [4"T] [2"9l [5] have indicated 
that the radio emission consists of core emission along 
the magnetic pole and one or more wider cones of emis- 
sion. The study of Arzoumanian et al.[3], as well as some 
earlier studies [ID], concluded that shorter period pulsars 
had more dominant core emission. Yet, several recent 
studies of the pulse polarization of young pulsars [TT] [28] 
and of pulsars showing three peaks in their profiles [19], 
suggest that pulsars with short periods, likely 7-ray pul- 
sars, have more dominant cone beams, implying a much 
wider beam of radio emission that in previous models. 

This paper will present the results of a study of both 
the slot gap and outer gap models, using the same popu- 
lation code and incorporating the revised radio emission 
model. We begin by reviewing the most widely studied 
models of pulsar high-energy emission, the polar cap 
and its extension to the slot gap, and the outer gap 
with its more recent extensions and modifications. We 
then discuss the traditional and revised models for ra- 
dio emission. The neutron star population synthesis is 
briefly discussed, since it is presented in more detail by 
Gonthier et al.[20] in these proceedings. The preliminary 
results of our study comparing high-altitude slot gap and 
outer gap models are then presented. Together with the 
study of the polar cap/low- altitude slot gap presented in 
Gonthier et al. 20J using the same set of evolved neutron 
stars, we can draw some important conclusions about the 
Geminga fraction in different models as well as the effect 



of radio and 7-ray emission geometry on the numbers of 
radio-loud and radio-quiet 7-ray pulsars. 



2 Gamma-Ray Emission Models 

Although the high-energy emission from rotation-powered 
pulsars has been studied for three decades, the particle 
acceleration and the location and mechanism of the emis- 
sion is still not understood. All models involve electro- 
static acceleration by an electric field parallel to the mag- 
netic field, but in different regions of the magnetosphere. 
All models also involve the production of electron-positron 
pairs and the cascades that determine the geometry of 
the observed 7-ray emission, by means of curvature and 
synchrotron radiation. 



2.1 Polar cap and slot gap 

In polar cap models, particle acceleration develops along 
open magnetic field lines above the neutron star surface. 
The accelerators divide into two types that depend on 
how charge is supplied and distributed. The two main 
subclasses are vacuum gap models [42] . where charges 
are trapped in the neutron star surface layers by bind- 
ing forces and a region of vacuum forms above the sur- 
face, and space-charge limited flow (SCLF) models [2\, 
where charges are freely emitted from the surface lay- 
ers. In SCLF accelerators, a voltage develops due to 
the small charge deficit between the real charge den- 
sity p and the Goldrcich- Julian charge density pcj — 
— 2e r2 • B (V • Ey = (p — pGj)ho), due to the curva- 
ture of the field [5] and to general relativistic inertial 
frame dragging [37] . where f2 and B are the rotation 
rate and surface magnetic field. Accelerated particles ra- 
diate 7-rays that create electron-positron pairs in the 
intense magnetic field and the En is screened above a 
pair formation front (PFF) by polarization of the pairs. 
The potential drop is thus self-adjusted to give parti- 
cle Lorentz factors around 10 7 , assuming a dipole mag- 
netic field. Above the PFF, force- free conditions could 
develop if the pair multiplicity is sufficient. This is likely 
for relatively young pulsars, which can produce pairs 
through curvature radiation [53], but older pulsars that 
can produce pairs only through inverse-Compton scat- 
tering are expected to be pair starved and their open 
magnetospheres would not achieve a force- free state [231 
[35] . Given the small radiation loss length scales for par- 
ticles of these energies, the high energy radiation will 
occur within several stellar radii of the surface. The ra- 
diation from electromagnetic cascades produces a hol- 
low cone of emission around the magnetic pole, with 
opening angle determined by the polar cap half-angle, 
PC (r) ~ (2-Kr/Pc) 1 / 2 , at the radius of emission r. 

More recent versions of the polar cap model[3_4J[35 
have explored acceleration in the 'slot gap' (Figure 1), a 
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Fig. 1 2D cross section of magnetosphere showing the loca- 
tion of the extended slot gap accelerator, for inclination angle 
a = 45°. 

narrow region bordering the last open field line in which 
the electric field is unscreened pQ. Near the open field 
line boundary, where the electric field vanishes, a larger 
distance is required for the electrons to accelerate to 
the Lorentz factor needed to radiate photons energetic 
enough to produce pairs. The PFF thus occurs at higher 
altitudes as the boundary is approached and curves up- 
ward, approaching infinity and becoming asymptotically 
parallel to the last open field line. If the electric field is 
effectively screened above the PFF, then a narrow slot 
surrounded by two conducting walls is formed. Pair cas- 
cades therefore do not take place near the neutron star 
surface in the slot gap, as do the pair cascades along 
field lines closer to the magnetic pole (core), but occur 
on the inner edge of the slot gap at altitudes of sev- 
eral stellar radii [34J. The high-energy emission beam of 
the slot gap cascade is a hollow cone with much larger 
opening angle than that of the polar cap cascade emis- 
sion. Even so, small values of both a and £, the angle of 
the magnetic axis and observer direction to the rotation 
axis respectively, (about 10°) are required to reproduce 
double-peaked profiles of observed 7-ray pulsars. 

Emission is also expected to occur from the primary 
particles that continue to accelerate in the slot gap at 
high altitudes because the potential in the slot gap is un- 



screened. Muslimov & Harding [35] have modeled the ac- 
celeration and emission pattern in the extended slot gap, 
finding that the electrons reach and maintain Lorentz 
factors ~ 10 . Such emission will include curvature, in- 
verse Compton and possibly synchrotron radiation. The 
emission pattern of the high-altitude slot gap radiation 
forms caustics along the trailing field lines and displays 
the geometry of the two-pole caustic model studied by 
Dyks & Rudak[l6]. The caustic pattern, as seen in Fig- 
ure[3J results from a near-cancellation of phase shifts due 
to aberration, timc-of-flight and curvature of the dipole 
field along the last open lines on the trailing side [33]. 



2.2 Outer gap 

Outer gap models :8„39] focus on regions in the outer 
magnetosphere that cannot fill with charge, since they lie 
along open field lines crossing the null surface, f2 ■ B = 0, 
where pcj reverses sign (see Figure [5]). Charges pulled 
from the polar cap, along field lines with little or no pair 
cascading from the polar cap accelerator, therefore can- 
not populate the region between the null surface and 
the light cylinder, and a vacuum gap forms. If outer 
gaps form, they can accelerate particles to high energy 
and the radiated 7 rays can produce pairs by interacting 
with thermal X-rays from the neutron star surface. Al- 
though the density of such X-ray photons is very small 
in the outer gaps, it is enough to initiate pair cascades 
since the newborn pairs accelerate in the gap, radiate, 
and produce more pairs. The gap size is limited by the 
pair cascades, which screen the gap electric field both 
along and across field lines, thus determining the emis- 
sion geometry. Young pulsars, having hotter polar caps 
and higher vacuum electric fields, tend to have narrow 
gaps stretching from near the null surface to near the 
light cylinder [7] while the gaps of older pulsars, having 
lower electric fields, are much thicker and grow with age 
[45] . When the gap fills the whole outer magnetosphere 
(at ages r > 10 7 yr) it ceases to operate, so that not all 
radio pulsars can emit 7 rays. Death lines in P-P space 
predict which pulsars can sustain outer gaps, depending 
on whether the X-ray photon field comes from cooling 
of the whole stellar surface or from polar caps heated by 
the energy deposited by the return flux of charges [46] . 

The recent outer gap model of Zhang et al.[46] has 
incorporated the dependence of gap geometry on a. The 
extent of the gap along the last open field line depends on 
a, because its inner radius, r in , which is assumed to be 
the null surface, is closer to the neutron star surface for 
oblique rotators. The growth of the gap across the field 
lines is limited by pair production on soft X-ray photons 
that originate from the neutron star due to surface cool- 
ing, polar cap heating or outer gap heating. Thus, the 
outer gap luminosity is expressed as 

Log = f 3 ((r),PB)L sd , (1) 
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Fig. 2 2D cross section of magnetosphere showing the lo- 
cation of the outer gap accelerator, for inclination angle 
a = 35°. 



where L sc i is the spindown luminosity, / is the size of the 
gap across field lines as a fraction of the total open field 
line volume, and is dependent on the average emission ra- 
dius in the gap, (r) , period P and surface magnetic field 
B. / is determined by the location of the pair formation 
front with respect to the last open field line, which is 
evaluated from the pair production condition 



ExE^l - cos9 Xl ) = 2(mc 2 f 



(2) 



where E x and E 1 are the energy of the soft X-ray and 
7-ray photon, and dx-y is the angle between their propa- 
gation directions along the field lines. We use this version 
of the outer gap model for the population synthesis pre- 
sented in this paper. 

Newer developments of the outer gap model include 
that of Takata et al. |43j . who solve Poisson's equation 
in two dimensions to determine the geometry of the gap 
both along and across field lines. They find that the gap 
may extend below the null surface in the presence of ex- 
ternal currents, as in the one-dimensional gap model of 
Hirotani et al. |26j . and that the gap width grows across 
field lines as the pulsar ages. This picture of the outer 
gap begins to resemble the slot gap in geometry, al- 
though the electrodynamics remains fundamentally dif- 
ferent. Cheng [6] has shown that this outer gap model 



may better reproduce light curves of known 7-ray pul- 
sars than the original model. 



3 Radio Emission Model 

Because the mechanism responsible for the radio beams 
is not understood, and more importantly because the ra- 
diation is coherent, it has not been possible to describe 
this emission using a physical model. The emission has 
therefore been described using empirical models, devel- 
oped over the years through detailed study of pulse mor- 
phology and polarization characteristics. The average- 
pulse profiles are quite stable and typically show a variety 
of shapes, ranging from a single peak to as many as five 
separate peaks. The emission is also highly polarized, and 
displays changes in polarization position angle across the 
profile that often matches the swing expected for a sweep 
across the open field lines near the magnetic poles in the 
Rotating Vector Model[38]. Rankin's[41j study of pulse 
morphology concluded that pulsar radio emission can 
be characterized as having a core beam centered on the 
magnetic axis and one or more hollow cone beams out- 
side of the core. Although Rankin's model assumes that 
emission fills the core and cone beams, other studies [31"] 
conclude that emission is patchy and only partially fills 
the core and cone beam patterns. 

Arzoumanian et al. (ACC)|3J fit average-pulse pro- 
files of a small collection of pulsars at 400 MHz to a 
core and single cone beam model based on the work of 
Rankin. The summed flux from the two components seen 
at angle to the magnetic field axis (modified by [18] to 
include frequency dependence v) is 



S(9, V) = -Fcoree 



where 
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and the index i refers to the core or cone, on is the spec- 
tral index of the total angle-integrated flux, Li is the 
luminosity of component and D is the distance to the 
pulsar. The width of the Gaussian describing the core 
beam is 



1.5°P- - 5 



(5) 



Pcore 

where P is the pulsar period in seconds. The annulus and 
width of the cone beam of Arzoumanian et al.[3], with 
frequency dependence of Mitra & Deshpande[32] are 
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The solid angles for the core and cone beams are 
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ACC found the ratio of core to cone peak flux to be 
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(10) 



where the frequency dependence has been added by [TB] • 
Thus in this model, F corc > F cono for nearly all pulsars 
and the core beam is completely dominant for short- 
period pulsars. 

We have incorporated several modifications to the 
ACC model for use in our population synthesis. The first 
change is to the ACC core-to-cone peak flux ratio result- 
ing from the recent work of Gonthier et al. [19| , who have 
carried out a study of 20 pulsars having three peaks in 
their average-pulse profiles, at three frequencies, 400, 600 
and 1400 MHz. They find a core-to-cone peak flux ratio 



25 P 

4 p-1.8 



1.3-0.9 
GHz ' 
-0.9 
'GHz ' 



P < 0.7s 
P > 0.7s 



(11) 



that is consistent with the core-to-cone peak flux ratio 
of ACC at periods above about 1 s, but predicts that 
pulsars with P < 0.05 s are cone dominated. Such a pic- 
ture is supported by polarization observations of young 
pulsars. Crawford et al.[TD][rT] measured polarization of 
a number of pulsars younger than 100 kyr, finding that 
they possess a high degree of linear polarization and very 
little circular polarization. Since conal emission typically 
shows high degrees of linear polarization, this strongly in- 
dicates that the emission from these young pulsars comes 
from part of a cone beam. Johnston & Weisberg[28]. 
studying polarization of 14 pulsars younger than 75 kyr, 
also find high degrees of linear polarization and flat posi- 
tion angle swings. They concluded that the emission was 
from a single wide cone beam, that core emission was 
weak or absent, and that the height of the cone emission 
is between 1% and 10% of the light cylinder radius. A 
high emission altitude for young (fast) pulsars was also 
found by Kijak & Gil[29], in their study of average-pulse 
profile widths and their dependence on frequency. As- 
suming that the edges of the pulse are near the last open 
field line, they find an emission radius of 
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where Tkg is in units of stellar radius. Our second modi- 
fication to the ACC model is therefore to add the radius 
dependence of cone beam emission from Eq. (Tl2|) above, 
so that Eq. ((6|) and J7J for the cone annulus and width 
are changed to 



9 = (1. - 2.63 8 w )p conc 

We — &w Pcone 

where 8 W = 0.18 [TJ5], and 



1.24° r^P 
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(13) 
(14) 

(15) 



The luminosities of the core and cone components are 
then 

j -^radio j -^radio 



1 + (r/r ) ' 



1 + (r /r) ' 



where 



^conc {pt-c 



l) l 



(16) 



(17) 



(18) 



^core ("cone + 1) T V50MHZ/ 

where a corc = —1.96 and a conc = —1.32, and 
L radio = 2.87 x 10 10 p- 13 P 0A mJy • kpc 2 ■ MHz 

as modified from ACC, where P is in units of 1 ss -1 . 

These changes to the radio emission model result in 
a significant revision of the geometry of radio beams in 
young and fast pulsars, from a narrow core beam emitted 
near the neutron star surface to a wide cone beam emit- 
ted at relatively high altitude. Such a revision produces 
greater visibility of the radio emission for fast pulsars and 
will have an influence on the predicted Geminga fraction. 



4 Pulsar Population Synthesis 

4.1 Galactic neutron star evolution 

The pulsar population synthesis code consists of two 
main parts. In the first part, neutron stars start from 
birth locations in the Galaxy, with given initial P, B, 
a, and space velocity, at constant birth rate. They are 
evolved through the Galactic potential and through their 
spin evolution to the present time. In the second part, 
each pulsar is assigned a radio and 7-ray flux and "de- 
tected" by different radio and 7-ray surveys. The first 
part of the code is described in some detail by Gonthier 
et al. |20) . In this paper, we have used the same set of 
evolved neutron stars to compute the numbers of radio- 
loud and radio-quiet pulsars expected assuming cither 
high-altitude slot gap or outer gap 7-ray flux and ra- 
dio flux predicted by the revised radio emission model 
described in Section ©. 

4.2 Emission geometry and luminosity 

To describe the emission geometry of each radiation model, 
we compute two-dimensional phase plots of emission in 
the (£, 4>) plane, where £ is the observer viewing an- 
gle measured from the rotation axis and <j> is the rota- 
tional phase. The method of computing the phase plots 
is identical to that of Dyks et al. .11] , who summed emis- 
sion tangent to field lines of a retarded vacuum dipole, 
taking into account relativistic effects of aberration and 
retardation. The emission was collected into bins on the 
sky (180 equal (j) bins, 180 equal ( bins), so that each 
2D element of the phase plot represents dL/d(dtfi, where 
dL = Idfl is differential luminosity and I and dfl are 
intensity and differential solid angle. Sample phase plots 
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Fig. 3 Plots of emission in the (£, (f>) plane for slot gap, outer gap and radio core and cone emission models, for P = 50 ms 
and a = 30° (top), 60° (middle) and 90° (bottom). 



for the slot gap, outer gap and radio cone and core beams 
are shown in Figure [3] Caustics on the trailing field lines 
are evident in the slot-gap and outer-gap phase plots. 
The shift of the cone beam relative to the core beam in 
radio phase plots is caused by aberration and retarda- 
tion, as has been observed in radio profiles [2"2"I[T5] . 



Slot gap and outer gap phase plots are computed 
for nine values a equally spaced between 0° and 90°. 
In our present calculation, the 7-ray phase plots depend 
only on a since we ignore any dependence of the radia- 
tion pattern on pulsar age. Radio phase plots are com- 
puted for the same values of a, but also for core and 
cone components separately and for five separate peri- 
ods P = 0.03,0.05,0.1,0.5,1.0 s and two different fre- 
quencies, v = 400 MHz and v = 1400 MHz. The mean 
intensity / in each phase plot is then normalized to the 
total luminosity 2L from both poles divided by the total 



emission solid angle 2fl in each model, where Q is the 
solid angle of emission from each pole: 



L 

n 



The phase-averaged flux, 



(vF v 



2nD 2 



(19) 



(20) 



is then determined by integrating the emission in the 
phase plot at a given £, after interpolation in a (and 
in P for the radio phase plots). This is equivalent to 
computing the phase-averaged flux in an observed pulse 
profile. 

For the slot gap model, we have assumed constant 
emission along the last open field line from the stellar 
surface to an altitude of 0.8 i?L, where i?L = c/(2 is the 
light cylinder radius. The total luminosity divided by 
solid angle from each pole is (from [34] ) 

LsG =e 7 [0.123 cos 2 a + 0.82 9 2 PC sin 2 a) ergs^sr" 1 



a 



SG 
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fsxlO 34 ( T5 ^ i75 )^P 5 / 7 , B<Q.1B C ^ 

X [2xl0 34 ( Tp ^ i75 ) 4/7 P 9 / 1 7 , B>0.1B cr 

where Pq.i = P/0.1 s and we have assumed an efficiency 
of conversion of primary particle energy to 7-ray emission 
of e 7 = 0.2. 

For the outer gap model, we have assumed constant 
emission along the field line having magnetic colatitudc 
£ = 9/6pc = 0.85, from the null surface to the light 
cylinder. The total luminosity of each pole is Log — 
/ 3 ((r) , P, B)L s d, as given above by Eq. and the solid 
angle of gap emission from each pole is ( [46] ) 




where b =< r > /i?isina. Using the above expressions 
we are able to reproduce the luminosity distribution in 
Fig. 2 of Zhang et al. [46]. 

For the radio emission model, all the field lines con- 
tained in the open volume are used and the flux S(6, v) 
from Eq ([3]) gives the differential luminosity in each sky 
bin: 

dV v = D 2 Si{9, v) sin Qd0dct> pc dv (23) 

emitted at altitude 1.8-R for the core component and at 
altitude given by Eq (fT2]) for the conal component, where 
4> pc is the magnetic azimuth. The total intensity of the 
phase plots are then normalized to 2 J J Si{6,v)dflD 2 . 



5 Predicted Populations of Gamma-Ray Pulsars 

Having assigned phase-averaged radio fluxes and 7-ray 
fluxes for either slot gap or outer gap models from the 
normalized phase plots to each evolved neutron star, we 
denote a source as radio loud if its radio flux exceeds 
the thresholds of any of the ten surveys described in 
Gonthier et al.[20j, and denote it as 7-ray loud if its 7- 
ray flux exceeds the flux thresholds of the 9 year EGRET 
or the GLAST 1 yr LAT (Large- Area Telescope) survey. 
The number of simulated neutron stars is normalized by 
matching the number of detected radio pulsars to the to- 
tal number detected by the ten surveys (which is equiv- 
alent to setting the neutron star birthrate). We use a 
revised EGRET threshold map, that takes into account 
the sky background model of Grenier et al.[2T], and a 
GLAST 1 yr LAT threshold map. In Table 1, we sum- 
marized the number of simulated radio-loud and radio- 
quiet 7-ray pulsars for slot gap and outer gap models 
predicted for EGRET and GLAST LAT telescopes, as 
well as the resulting Geminga fraction in parentheses. In 
both slot gap and outer gap models, the number of radio- 
quiet 7-ray pulsars far exceeds the number of radio-loud 
7-ray pulsars, thus predicting a large Geminga fraction. 
Models where the high-energy emission occurs mostly in 
the mid- or outer magnetosphere therefore produce many 



Geminga-like pulsars and very few radio- loud pulsars. 
The Geminga fraction in the outer gap model is signif- 
icantly higher than in the slot gap model, with too few 
radio-loud pulsars to account for the number EGRET 
detected and only 1 predicted to be detectable in 1 year 
by the GLAST LAT, even though the total number of 
detectable 7-ray pulsars is a factor of 3 higher. The rea- 
son for the much higher number of 7-ray pulsars in the 
outer gap model is clear from Fig. [4] which shows the flux 
distribution of the pulsars detectable by both models. It 
is evident that there are many more pulsars with high 
7-ray fluxes in the outer gap model, although the dis- 
tributions peak at about the same flux in both models. 
This is because although the luminosity distributions are 
similar, the outer gap pulsars have smaller effective solid 
angles, as shown in Figure [5] The age distributions of 
the detected EGRET and LAT 7-ray pulsars, as shown 
in Figure [6] are quite similar for the slot gap and outer 
gap and the distribution strongly peaks at small ages. 
However, the radio-loud pulsars appear to be somewhat 
more spread-out in age, although the statistics are lim- 
ited. The Galactic latitude distribution in Figure [7] shows 
that both the slot gap and outer gap 7-ray pulsars are 
concentrated near the Galactic plane, with the distri- 
bution of outer gap pulsars being much more strongly 
peaked at \b\ = 0. The |6| distribution of slot gap pulsars 
more closely resembles the latitude distribution of the 
EGRET sources. 

Gonthier et al. [20] have presented results for the num- 
bers of radio-loud and radio-quiet pulsars with detectable 
emission from the low-altitude slot gap cascades, as well 
as from high-altitude slot gap emission. Those results 
are thus complimentary to these, and the results for 
the high-altitude slot gap can be compared. Both the 
high-altitude and low-altitude emission is expected to be 
present for pulsars having slot gaps, so the numbers in 
Table 1 should be combined to give the total number of 
7-ray pulsars expected for the polar cap/slot gap model. 
The Geminga fraction for the low-altitude slot gap (0.28 
for EGRET and 0.35 for GLAST) is much lower than 
for the high-altitude slot gap. This is because the low- 
altitude slot gap cascade emission occurs at 3-4 stellar 
radii, much closer to the radio emission altitude for most 
of the 7-ray loud pulsars, especially the ones that are 
core dominated. The low-altitude slot gap also produces 
enough EGRET radio-loud pulsars and can account for 
some fraction of the unidentified 7-ray sources as radio- 
loud 7-ray pulsars, whereas we found that the outer mag- 
netosphere models do not produce enough EGRET pul- 
sars. 

Jiang & Zhang [27] have studied statistical properties 
of 7-ray pulsars in the outer gap model. From their calcu- 
lation, they predict that 8 radio-loud and 24 radio-quiet 
7-ray pulsars are detectable by EGRET, and 78 radio- 
loud and 740 radio-quiet 7-ray pulsars are detectable by 
GLAST. The predicted Geminga fractions are therefore 
0.75 and 0.9 for EGRET and GLAST. These fractions 
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Fig. 4 Flux distribution of detected 7-ray pulsars for the slot gap and outer gap models for EGRET (left) and LAT (right). 
Radio-loud pulsars are in dark shading and radio-quiet pulsars are in light shading. 
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Fig. 5 log(L 7 /i ? 7 D 2 ) distribution of detected 7-ray pulsars for the slot gap and outer gap models for EGRET (left) and 
LAT (right). 
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Fig. 6 log(age) distribution of detected 7-ray pulsars for the slot gap and outer gap models for EGRET (left) and LAT 
(right). 
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Fig. 7 Galactic latitude distribution of detected 7-ray pulsars for the slot gap and outer gap models for EGRET (left) and 
LAT (right). 
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Table 1 Predicted numbers of 7-ray pulsars and Geminga fraction" 

Slot gap Outer gap 

Low altitude 6 High altitude 

Instrument Radio-Loud Radio-Quiet Radio-Loud Radio- Quiet Radio-Loud Radio-Quiet 
EGRET 25 10 [028) 3 31 JOM) 0^1 118 (099)" 

GLAST LAT 94 51 (0.35) 13 128 (0.91) 1 258 (0.99) 

"in parentheses 

6 From Gonthier et al.[20] 



are significantly lower than our calculated Geminga frac- 
tions for the outer gap. But there are at least two ma- 
jor differences in our calculations. First, we have used 
more recent full-sky threshold maps for EGRET and 
the GLAST LAT, taking into account the sky coverage 
of each telescope and detection above the intense inter- 
stellar background, which has the effect of reducing the 
number of detected 7-ray pulsars. Jiang & Zhang [27) 
used only in-plane and out-of-plane sensitivities, which 
do not accurately account for the large variation in sensi- 
tivity with Galactic coordinates. Second, we have taken 
into account the relative spatial orientations of the radio 
and 7-ray beams, as well as their solid angles, whereas 
Jiang & Zhang have taken into account only the different 
solid angles and beaming fractions of the radio and 7-ray 
emission. This is important since the outer gap emission 
direction is generally at large angles to the radio emis- 
sion direction, and there is a large region of phase space 
where radio and 7-ray beams are not visible to the same 
observer, e.g. for a<40°. 



6 Conclusions 

We have determined the numbers of radio-loud and radio- 
quiet 7-ray pulsars using two different 7-ray models, from 
the same set of evolved Galactic neutron stars and us- 
ing the same radio emission model. The radio emission 
model takes into account the most recent studies of radio 
pulse emission morphology and polarization, which gives 
wide, relatively high altitude radio cone beams for the 
young 7-ray bright pulsars. The full geometry of both ra- 
dio and 7-ray emission is modeled, including relativistic 
effects of retardation and aberration, and distortion of 
the open field lines due to a retarded dipole description. 
The Geminga fraction is large for models such as the 
extended slot gap and outer gap where 7-ray emission 
occurs at high altitude in the pulsar magnctosphere, for 
a radio beam model that describes the bulk of the ra- 
dio pulsar population. Even the larger radio beams of 
young pulsars (P < 50 ms) emitted at high altitude, pro- 
duce few radio-loud pulsars and a low Geminga fraction, 
since the number of fast pulsars (P < 50 ms) with very 
large radio beams is not large enough. For the bulk of 
the pulsars with P > 50 ms, the radio beam size decreases 
rapidly while the 7-ray beams remain large. Both slot 
gap and outer gap models yield a large spread of the 



/ 7 -D 2 vs. Lsd distribution. In general, the radio- loud 
pulsars are closer and have larger L^/F^D 2 than the 
radio-quiet pulsars. From the results of this paper for 
high-altitude 7-ray emission models, combined with the 
results of Gonthier et al.[20] in this proceedings for the 
polar cap/low-altitude slot gap model, we can come to 
the important conclusion that if many of the EGRET 
unidentified sources are radio-loud 7-ray pulsars (and 
not pulsar wind nebulae), the 7-ray emission must come 
from relatively low altitudes. Of course, it is quite possi- 
ble that different emission models may apply to pulsars 
of different ages. Future observations with GLAST of the 
number of radio- loud 7-ray pulsars, and a related limit 
on the Geminga fraction, will be able to distinguish be- 
tween low- altitude and high-altitude emission models. 
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